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Abstract: A nonclassical light source is used to demonstrate experi- 
mentally the absolute efficiency calibration of a photon-number-resolving 
detector. The photon-pair detector calibration method developed by Klyshko 
for single-photon detectors is generalized to take advantage of the higher 
dynamic range and additional information provided by photon-number- 
resolving detectors. This enables the use of brighter twin-beam sources 
including amplified pulse pumped sources, which increases the relevant 
signal and provides measurement redundancy, making the calibration more 
robust. 
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1. Introduction 

Quantum optics enables one to make measurements that are more precise than the fundamental 
limits of classical optics H]. Central to this capability are quantum optical detectors, those that 
are sufficiently sensitive to discern the inherent discreteness of light. These detectors are key to 
emerging quantum technologies such as quantum imaging and lithography, in which the stan- 
dard wavelength limit to resolution is surpassed by using quantum states of Ught and photon- 
number sensitivity. However, the majority of quantum optical detectors have a response that sat- 



urates at only one photon, imposing a significant limitation on the brightness of the optical fields 
that can be used for such quantum technologies. The result of this binary detector response is 
a measurement that can discriminate only between zero photons and one or more photons ar- 
riving simultaneously; the detector produces an identical response for any number of photons 
greater than zero. To allow increased brightness of the light sources used in the technologies 
outlined above (and the concomitant improvement in accuracy that this brings), it is necessary 
to use detectors that can discern the number of photons incident simultaneously on the detector 
— a photon-number-resolving detector (PNRD). Indeed, the development of PNRDs is an ac- 
tive area of research including photomultiplier tubes, the extension of avalanche photodiodes 
(APDs) to higher photon number through time multiplexing G] 111 or spatially-multiplexed ar- 
rays m, transition-edge sensors (TES) |l5]|6l|7l, charge-integration photon detectors (CIPS) El, 
superconducting single-photon detectors |I9][10|, visible-light photon counters (VLPCs) ifTTI . 
and quantum dot field effect transistors lfT2l[T3l . 

Calibration of optical detectors is a difficult problem. The standard approach uses a previ- 
ously calibrated light source. The drawback of this approach is that errors in the source bright- 
ness translate directly into errors in the detector efficiency calibration. The converse is also 
true, leading to a detector and light-source calibration dilemma. To get around this, brightness 
calibration is typically based on a fundamental physical process, for example, the luminosity of 
blackbody radiation of gold at its melting point ll4!j, or heating of a cryogenic bolometer ifTsl . 
Such methods are suitable for calibrating bright light sources. In contrast, using this method 
to calibrate detectors operating at the quantum level (i.e. fields containing only a few photons) 
requires sources with powers on the order of a femtowatt to avoid saturation. Such sources are 
impractical. 

Nonclassical states of light allow us to circumvent the horns of the dilemma due to their 
behavior in the presence of loss. Realistic detectors can be modeled by optical loss (i.e. at- 
tenuation) followed by unit efficiency detectors. The transformation of classical states (e.g. a 
coherent state, thermal state, etc.) under loss is parameterized only by source brightness, which 
scales linearly with the loss. In contrast, nonclassical states change their statistical character 
upon experiencing loss: correlations in optical phase, intensity, photon number, and electric 
field transform in a non-trivial way under loss. Based on this fact, Klyshko proposed a way to 
calibrate detectors based on the statistical character of light rather than its brightness lfT6l . This 
approach relies on spontaneous parametric downconversion (SPDC) as a light source in which a 
photon is simultaneously created in each of two optical modes, usually denoted signal and idler 
Since the photons are created in pairs the two output modes are perfectly correlated in photon 
number. Thus detection of the idler without the simultaneous detection of the signal photon can 
only be caused by loss in the signal arm. Measuring the detected rate of idler photons /?, and 
photon pairs Rc then allows a calibration of the detector efficiency Tjj 



and vice versa with ; ^ s. This efficiency estimation was shown experimentally in lfT7lll8l[T9l 

The Klyshko scheme is limited by three factors. First, it relies on the implicit assumption 
that at most one photon pair is emitted at a time, a feature which can be violated if the SPDC 
is pumped strongly. Indeed more than one pair is often desirable, as in continuous-variable 
experiments. In this case, the detector efficiency can be obtained by first lowering the pump 
power and using the Klyshko method 1211 . However, the detector is now calibrated outside the 
regime of its intended use, which could necessitate subsequent assumptions such as the inde- 
pendence of the estimated efficiency from the pump power. Thus, the direct in situ calibration 
of detectors would be desirable. Second, the Klyshko scheme is primarily designed for single- 



photon detectors and is not directly translated to PNRDs. And third, this calibration method is 
highly sensitive to the input state quality and measurement uncertainties. Because the number 
of measurements taken is exactly the number needed to determine the efficiency, any errors in 
the measurements will propagate directly into the efficiency estimation. 

Here we present a technique for measuring the absolute quantum efficiency of a detector 
based on the Klyshko method, but explicitly taking into account multiple photon events. This 
improves the calibration accuracy and allows for calibration of a PNRD. This approach utilizes 
the PNRD capability to measure the photon-number distribution of an optical mode. Using two 
PNRDs the joint photon-number statistics between the two electromagnetic field modes, in- 
cluding photon-number correlations and individual photon-number distributions, of the SPDC 
source can be determined. For each element of the resulting joint photon statistics, one can 
find a formula giving the detector efficiencies of the two PNRDs. For the zero- and one-click 
elements, this reproduces the Klyshko result. However, the increased dynamic range of PNRDs 
allows the use of brighter sources that produce measurable rates in the higher photon-number 
elements of the joint statistics. We then use optimization techniques to estimate the detector 
efficiencies from the increased number of measurements. This added redundancy improves the 
tolerance of our scheme to background light and statistical noise. We experimentally demon- 
strate this efficiency estimation method with two time-multiplexed PNRDs [Jl. 

We begin by introducing a general treatment of PNRDs and then use this as a basis for de- 
scribing our generalized Klyshko method. This is followed by the description of an experiment 
to test the efficiency estimation with PNRDs. 

2. Photon-number-resolving detectors 

Photon-number-resolving detectors (PNRDs) are a class of photodetectors that have a unique 
response for every input photon-number state within their range. Ideally these responses can be 
perfectly discriminated. However, the less than perfect efficiency of realistic detectors causes 
these responses to overlap, and thus does not allow for direct photon-number discrimination. 
Overlap of detector responses can also arise from the detector electronics (e.g. amplification) or 
the underlying detector design. Despite this overlap, the linear relationship between the detector 
response and the input state allows for the reconstruction of the input photon statistics from the 
measured outcome statistics. This linear relationship is encapsulated by 

P„ = Tr [pn„] , (2) 

where p is the input-state density matrix, P„ is the probability for the nth measurement outcome 
and n„ is the associated positive operator-value measurement (POVM) operator. Consider the 
matrices expressed in the photon-number basis. Since PNRDs do not contain an optical phase 
reference, the off-diagonal elements of n„ are zero, meaning the photon-number-resolving de- 
tection is insensitive to off-diagonal elements in p . It is thus useful to write the diagonal ele- 
ments of p, the photon-number statistics, as a vector a. Similarly, we write the outcome prob- 
abilities {Pi} as a vector p. In the following, we truncate a at photon number N ~\, where 
is the number of detector outcomes, although this is not strictly necessary. 

The POVM operators of a general PNRD can be modeled by dividing the detector perfor- 
mance into two components: efficiency and detector design, described by matrices F and L 
respectively. In the photon number basis, a POVM element can be written as fl24l 

n„=^[F.L(77)]„,Jm)(m|, (3) 

m 

where [F • L ,„ corresponds to the probabilty of detecting « out of m incidence photons. As 
pointed out in the introduction, detector efficiency Tj can be modeled by a preceding optical loss 



/ of (1 — Tj). In the context of a PNRD, loss causes the photon-number statistics to transform 
according to a ^ L (77) a, where 



These matrix elements transform the state by lowering photon numbers from j to /, representing 
the loss of photons through a binomial process with probability / = (1 — 77). 

Although the detector-design component of the model depends on the detailed functioning 
of the device, a large class of PNRDs - mode-multiplexers - can be treated in the same way. 
These detectors divide an input optical field mode into many spatial and/or temporal modes 
and then use single-photon detection on each mode to achieve number resolution. Examples 
include nanowire superconducting detectors ||9l, VLPCs ifTTl . intensified charge-coupled de- 
vices (CCDs) E2I, integrated APD aiTays ElllSl, and time-multiplexed detectors (TMD). All 
of these detectors suffer from detector saturation; the one-photon detector response occurs if 
two or more photons occupy the same mode. This saturation effect is modelled by a detector 
design matrix F = C, the "convolution" matrix, which has a general but complicated analytic 
form given in 1241 (in the context of the TMD). The form of C depends on relatively few pa- 
rameters comprising the splitting ratios of the input mode into each of the multiplexed modes 
and the total number of these modes. 

As an example, in a CCD array detector the pixel shape and size defines the detected optical 
mode. The spatial overlap of the detector mode of each pixel with the incoming optical mode 
then gives the corresponding splitting ratio. The total number of pixels is the total number of 
multiplexed modes. As the number of multiplexed modes goes to infinity the C matrix goes to 
the identity. All PNR detectors can be described by a POVM, which can be reconstructed by 
detector tomography 12411251 . For PNR detectors that do not rely on mode multiplexing, such 
as transition-edge superconducting detectors, and single APD detectors, one must factor the 
POVM elements into an F matrix and an L matrix to apply the calibration procedure described 
in this paper. 

Focusing on the specific case of mode-multiplexed detectors, one can rewrite the Eq. ^ as. 



where the elements of the rth row of C • L (77) are the diagonals of the ;th operator in the POVM 
set {ft,}. 

With PNRDs in two beams, denoted 1 and 2, one can measure not only the individual photon- 
number statistics di and 02, but also the joint photon-number distribution of these two beams. 
This distribution is written as the joint photon statistics matrix C7, where a,„,n is the probability 
of simultaneously having m photons in mode 1 and n photons in mode 2. We extend Eq. (|5]l to 
relate the probability P,„,,„ of getting outcome m at detector 1 and outcome n at detector 2, to 
the joint photon statistics a, 



where subscripts indicate the relevant detector and is the transpose of A. Joint photon statis- 
tics are a measure of photon-number correlations in beams 1 and 2, and are thus sensitive to loss 
as discussed in the introduction. We use this description of PNRDs to generalize the Klyshko 
calibration. 

3. Generalizing the Klyshko method 




p = C-L{ri)-a, 



(5) 



P = Ci-L(t7i).(J-LT(772)-CJ, 



(6) 



The assumption that only a single photon pair is generated at a time in the Klyshko efficiency- 
estimation scheme is only valid for very low SPDC pump powers. We can determine in what 



manner the Klyshko scheme breaks down when more than one pair is produced. Simply assum- 
ing an additional contribution of pairs (72,2 to the pair rate (7i i changes the efficiency estimate 
of Eq. ([B to 

R, (7i.i + (77,T7,-277,-277,+4)(T2.2 

f?i = 5- = -7^ ^ fJi- (7) 

This overestimates the efficiency for (72,2 7^ 0. Similarly, sensitivity to the detailed photon statis- 
tics of the input states acts to degrade the accuracy of the efficiency estimate. Let us consider 
the effect of background light (e.g. fluorescence from the optical elements and detector dark 
counts) on the efficiency estimation of Eq. ([T]i. Background photons are uncorrected between 
the signal and idler beams, increasing the singles rate, and thus increasing /?,-. This will make 
the estimated efficiency of the detector fjs lower than the actual efficiency rj^. 

To generalize the Klyshko scheme to PNRDs and avoid the above limitations we need to 
determine the true state generated by SPDC photon-pair sources. Ideally, these sources produce 
a 'two-mode vacuum squeezed state' 



where A is proportional to the pump beam energy, and \n)^(i) is an n-photon state of the signal 
(idler) mode. Having only one free parameter, it is tempting to use this state in the generalized 
efficiency estimation. However, this state can only be generated with careful source design 1261 . 
Instead, photons are typically generated in many spectral and spatial modes in the signal and 
idler beams |[27l . Depending on the number of modes in the beams, the thermal photon-number 
distribution in Eq. ^ changes continuously to a Poisson distribution 1281. Consequently, it 
would be incorrect to assume the source produces an ideal two-mode vacuum squeezed state. 
Still, the number of photons remains perfectly correlated between the two beams. Without 
access to the number of generated modes we can make only the following assumption about 
the joint photon statistics of the source 

^in.n ^ ' ^m^rl: (9) 

where {c, } are arbitrary up to a normalization constant and 5,„,„ is the Kronecker delta. 

Inserting the joint photon statistics defined by Eq. (|9]l into Eq. (|6]l we arrive at the basis for our 
generalized Klyshko method. Since Ci and Ci of the detectors are known the predicted joint 
outcome probabilities P are highly constrained having A^^ elements uniquely defined by the 
parameters in {c,} and the two efficiencies rji and 772. Consequently, a measurement of the 
joint outcome statistics specifies Tji and TJ2 with a large amount of redundancy; the efficiencies 
are overdetermined. In order to correctly incorporate all measured outcome statistics into the 
efficiency estimates a numerical optimization approach is used. We minimize the difference G 
between the measured outcome statistics R and the predicted outcome statistics P, which are 
determined by {c, } , Tji, and 772 

G = R-Ci-L(77i).(J.lT(772)-CJ. (10) 

r r 21) '/^ 

This is done by minimizing the Frobenius norm F = <Tr (G) > to find the optimal 771 
and 772 - our estimates of the PNRD efficiencies. Using the Frobenius norm makes this a least- 
squares optimization problem over {c, } , 771 and 772, where < 77, < 1. However, this method of 
efficiency estimation is amenable to other optimization techniques such as maximum-entropy 
or maximum-hkelihood estimation. We use the Matlab® function Tsqnonneg' (with the con- 
straints G„,j„ > 0). 




r|i, Detector efficiency 

Fig. 1. A typical optimization residual (F) for simulated joint outcome statistics as a func- 
tion of the two PNRD efficiencies T]i and rjo. There is only one minimum, suggesting that 
the problem is convex. 



This efficiency estimation is similar to the optimizations performed in state or process tomog- 
raphy, which are known to be convex problems (i.e. there are no local minima) 1291 . However, 
we are now estimating parameters in both our state p and the POVM set {ft/} . Because this 
is not guaranteed to be convex ||29ll , we simulated a variety of measured statistics to test for a 
single minimum. Using Ci and C2 for the PNRDs in our experiment (TMDs), with several sets 
of photon statistics {c, } , and a range of efficiencies 771 and 772, we simulated various outcome 
statistics. In all cases, the optimization reproduced the correct efficiencies and only a single 
minimum was observed. Figure [T] shows the result of a typical simulation displaying a global 
minimum in the optimization residual (the minimum value of F for a given pair of efficiencies 
rji and T72). This is a good indication that efficiency estimation is a convex problem. 

4. Experimental setup 

To experimentally demonstrate and test our efficiency estimation method, time-multiplexed de- 
tectors (one possible realization of a PNRD) were employed. In a TMD, the input optical state 
is contained in a pulsed wavepacket mode. The pulse is split into two spatial and several tempo- 
ral modes by a network of fiber beam splitters and then registered using two APDs [J]. APDs 
produce largely the same response for one or more incident photons. The TMD overcomes this 
binary response by making it likely that photons in the input pulse separate into distinct modes 
and are thus individually registered by the APDs 13] . The TMD is a well-developed technol- 
ogy, which makes this an ideal detector to test our approach to detector efficiency estimation. 
The convolution matrix C for this detection scheme is calculated from a classical model of the 
detector using the fiber splitting ratios |3i, and is also reconstructed using detector tomography 




Fig. 2. The experimental setup. A KDP crystal two-beam source is pumped by an amplified 
Ti:Sapph laser. The two generated beams propagate collinearly and are orthogonally po- 
larized. Each beam is measured by a time-multiplexed photon number resolving detector. 
Details are given in the text. 



Il24l |251 . Loss effects in TMDs have also been thoroughly investigated ll30ll . In our experi- 
ments, the TMDs have four time bins in each of two spatial modes, giving resolution of up to 
eight photons, with a possible input pulse repetition rate of up to 1 MHz. Field-Programmable- 
Gate-Array (FPGA) electronics are used to time gate the APD signals with a window of 4 ns, 
which significantly cuts background rates. The joint count statistics R are accumulated by the 
electronics and transferred to a computer for data analysis. 

The experimental setup is centered on a nearly-two-mode SPDC source 1261 as depicted 
in Fig. [21 The twin beam state produced by SPDC in a potassium dihydrogen phosphate 
(KDP) nonlinear crystal consists of two collinear beams with orthogonal polarizations. The 
pulsed pump (415 nm central wavelength) driving the SPDC is a frequency-doubled amplified 
Ti:Sapphire laser operating with a 250 kHz repetition rate. A pick-off beam is sent to a fast 
photodiode (PD) that is used to trigger the detection electronics. Dichroic mirrors (DM) and a 
red-pass color glass filter (RF) are used to separate the blue pump from the near-infrared (830 
nm central wavelength) SPDC light. A polarizing beam splitter (PBS) is used to separate and 
direct the two co-propagating downconversion beams into separate single-mode fibers (SMFs) 
connected to the time-multiplexed detectors (TMDl and TMD2). The joint statistics R of the 
two TMDs are recorded for a range of pump powers between 1 and 55 mW in order to estimate 
the two TMD efficiencies at each power. To examine the spectral response of the detector effi- 
ciency, one could tune the wavelength of the SPDC source by adjusting the pump wavelength 
and the crystal orientation. 
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Fig. 3. The estimated detector efficiencies for TMDi (■) and TMD2 (A) as a function of 
tlie average SHG power pumping tlie SPDC. Three distinct regimes (regions 1, 2, and 3) 
are indicated by shading. 



5. Estimated efficiencies 

For each pump power we determine the optimum detector efficiencies that are consistent with 
the measured statistics at both PNRDs, shown in Fig. [3] Three different regimes are observed. 
At low powers (up to 6 mW) the estimated efficiency increases with power Between powers of 
6 and 40 mW the estimates appear constant. At 40 mW there is a sudden jump in the estimated 
efficiency (approximately twice the previous value); above this power the estimates remain 
constant. Continued investigation revealed that the second-harmonic generation (SHG) pro- 
cess qualitatively changed its behavior at 40 mW: the increased pump power induced unwanted 
higher-order nonlinear effects, resulting in the generation of additional frequency components 
other than the second harmonic and a change in the spatial mode structure. This changed both 
the transmission of the short wave pass filter (DM and BF) and the efficiency of the fiber cou- 
pling into our detectors. Since this behavior is outside the scope of our investigation we omit 
this data from further discussion. 

The TMD efficiency should be independent of the average photon number of the state and 
thus the pump power. This is true in the second region of Fig.[3]but not the first. By reconstruct- 
ing the joint photon-number distribution of the input state (<7,„,„) using the estimated efficien- 
cies, one can gain insight into the estimation accuracy of the detector efficiency. This serves 
as a partial check for our assumption that the number of photons in the two beams is equal. 
In Fig.m we show the reconstructed joint photon statistics for two pump powers in regimes 1 
and 2. In the second regime, the photon-number distribution is largely diagonal: only 10% of 
the incident photons arrive without a partner in the other beam. In contrast, the state in the low 
power regime has significant off-diagonal components, with 43% of the photons arriving alone. 
This suggests that at low powers the reference state is corrupted by background photons, possi- 




Fig. 4. The reconstructed photon statistics o for m and 112 photons in beam 1 and 2, respec- 
tively, measured at a pump power of (a) 1.5 mW and (b) 30 mW. The presence of significant 
off-diagonal elements indicate that the input state is corrupted by background. 



bly fluorescence from optics in the pump beam path, pump photons leaking through our filters, 
or scattered pump photons penetrating our fiber coatings. Contributions from dark counts are 
expected to be negligible, since the specified dark count rates of our detectors are significantly 
lower than these other effects. In the next section we attempt to remove this background in 
order to better estimate the efficiency. 

6. Compensation for background light 

We investigate two methods of dealing with background light. In the first, a parameterized 
background contribution is incorporated into the efficiency estimation procedure. The second 
attempts to measure the outcome distributions due to the background alone and then subtract 
these from the efficiency data. 

As pointed out in Section 3, the efficiency estimation is greatly overdetermined. One ex- 
pects that a modeled background could be added to our input state model, Eq. without 
jeopardizing the convergence of the optimization. This background would be entirely uncor- 
rected, possibly making its contribution to the joint outcome statistics easily distinguishable 
from the SPDC contribution. We model the photon-number statistics of the background in each 
beam by a Poisson photon-number distribution d{n) = a"exp(— a)/«! Ol], which is fixed by 
a single parameter - the average background photon number, a = (ng). Consequently, only 
two additional parameters (one for each beam) enter into the efficiency estimation, keeping 
it overdetermined. Unfortunately, we theoretically found that the loss in one beam transforms 
the outcome statistics in a manner similar to background in the other beam. Thus, the prob- 
lem is no longer convex; there is a set of equally optimal points {((«b) , rj)}. To show this we 
compare two different two-mode number-coiTelated states cj'*'^', similar to the state in Eq. 
that undergo the addition of uncorrected background light and loss respectively. The addition 
of background to beam 1 of state a'* is given by the convolution of the Poisson distribution 
with input state, (defined by arbitrary {c^}). The elements of the background-added joint 
photon-number statistics a'* are given by the sum of all possible ways to add photons from 
the Poisson background, d, to the first beam of the initial state, ff'*, that add up to a particular 



number of photons in the final state, 



k 

L^H-^-w- (11) 

m=0 

This can be written in terms of a matrix product 

((««)) •(T^ (12) 

where the elements of D,„,„ = c/(n — m) are the probabilities of having an additional 

n — m photons from the Poissonian background. To compare this background-added case with 
the situation in which there is no background, but there is loss, we assume a loss / = (1 — T]) 
in beam 2 of the second state (7^ (defined by arbitrary {c^})- We then attempt to show that for 
some and (each with perfectly correlated photon statistics) the two resulting states are 
equal, that is 

Z)((«fi)).(j'* = CJ«.LT(77). (13) 

EUmination of {c^} and {c^} from the resulting equations is facihtated using a computer al- 
gebra program, and we find that there is a range of / and (ng) < 1 that solve Eq. ( fTSl ), indicating 
that it is not possible to simultaneously fit for background and efficiency. This emphasizes the 
fact that one cannot distinguish between loss and Poisson background for the number-correlated 
states. Figure |5] shows / as a function of (ng) for PNRDs with a maximum photon number of 
M = 1 to 20. The standard Klyshko case corresponds to M = 1 . Note that as M becomes large 
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Fig. 5. The loss / in beam 1 of a twin-beam state ct* that results in the same joint outcome 
statistics as the addition of a Poissonian background, with average photon number (ns) , 
to beam 2 of another twin-beam state, o'^, for some and CT*. M is the photon number 
range of PNRDs in beams 1 and 2. 



the loss curve converges to line of slope 0.068 suggesting that even PNRDs with an infinite 
photon number range would not allow one to distinguish loss and background. This suggests 
that another approach to background light should be used. 

Another approach to address the background contribution to the efficiency estimation at- 
tempts to subtract an independently measured background contribution from the click statis- 
tics. For each pump power the pump polarization is rotated by 90 degrees, extinguishing the 
SPDC and allowing the measurement of the joint outcome statistics due to background alone. 
Generally, the statistics of two concurrent but independent processes is the convolution of the 
statistics of the two processes. However, the measured outcome probabilities Pm of both light 
sources combined is not a simple convolution of the background outcome probabilities Pb and 
the twin-beam outcome probabilities P5 

Pm^Ps*Pb. (14) 




Fig. 6. The reconstructed photon statistics a for «i and 772 photons in beam 1 and 2, re- 
spectively, after subtracting an independently measured background. For pump powers (a) 
1.5 mW and (b) 30 mW, as in Fig. |4] The reduction of significant off-diagonal elements 
indicates that the background subtraction method works. 



This is due to the fact that there is an effective interaction between the background and 
twin-beam signal detection events due to the strong detection nonlinearity of the APDs (i.e. 
the detectors saturate at one photon). However this saturation effect can be eliminated by first 
applying the inverse of the C matrices to the measured statistics 

p;, = c^'PM(c|)-^ (15) 

P^ = C^'Ps(C|)-^ (16) 

which gives the estimated photon-number statistics prior to the mode multiplexing. These non- 
mode-multiplexed statistics can be assumed independent so that, 

p;, = [crip,(cir']*p^. 

We then use the convolution theorem to find 



where ^ indicates the Fourier transform, and the matrix division is element by element. Using 
Ps we estimate the efficiency as previously by the method in Section 3. 

To test the accuracy of this background subtraction method we reconstruct the joint statistics 
at the same powers as in Fig. [3] We find that in both cases the off-diagonal components are 
significantly reduced. In the low power regime, now only 16% of incident photons are not 
part of a pair, and at higher power this becomes just 4%. With the background subtracted, the 
estimated efficiencies are plotted in Fig. |7] and are now in better agreement with the expected 
constant detector efficiency through the first two regions. However, the estimates still drop off 
as the power goes very low. Although we measure the background, we do not do so in situ; we 
need to change the apparatus by rotating a half wave plate (HWP). Thus, we are not guaranteed 
that this is equal to the background present during the efficiency estimation. Furthermore, errors 
in background measurements are more significant at low powers as background then forms a 
larger component of the outcome statistics. 

Also plotted in Fig.|2]is the Klyshko efficiency. In contrast to the increased dynamic range of 
our method the standard Klyshko efficiency increases with pump power, evidence that higher 
photon numbers in the input beams distort the estimated efficiency. 

The average efficiency across the second region was found to be 9.4% ± 0.4% for detector 1 
and 8.0% ± 0.4% for detector 2, where the errors are the standard deviations. These relatively 




10 20 30 

Pump power (mW) 



Fig. 7. The estimated detector efficiencies for TMDi (■) and TMD2 (A) determined from 
background subtracted outcome statistics, plotted as a function of the average pump power. 
Also plotted is the Klyshko efficiencies that would have been estimated for single photon 
detector 1 ((}) and 2 (Q)- The standard Klyshko method overestimates the efficiencies for 
high powers. The dotted lines indicate the average efficiencies of the two PNRDs. 



low efficiencies are only partly due to the quantum efficiency of the avalanche photodiodes 
themselves, which is specified to be 60% ±5% at our wavelength. Bulk crystal SPDC sources 
ordinarily emit into many spatial modes, which makes coupling into a single-mode fiber diffi- 
cult and inefficient. Typical coupling efficiencies are less than 30% 1321. 
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Fig. 8. The reconstructed average photon number (+) as a function of the pump power. 
SPDC theory predicts a linear relationship. The line is a linear fit. 



As a final check of the reconstructed photon statistics, we plotted the average reconstructed 
photon number as a function of pump power in Fig. [8] The relationship is linear, as expected 
when taking into account the higher dynamic range of a TMD in comparison to standard APDs 

m- 

7. Conclusion 

Despite being one of the seven base SI physical quantities, the working standards for luminous 
intensity have a relative accuracy of only 0.5% 1341 , compared to 10^^^ and better for the sec- 
ond ||35l . Relying on a light beam of a known intensity, the efficiency calibration of detectors 
is similarly limited. Quantum states of light give us the opportunity to bypass the working stan- 
dards and calibrate detectors directly. Indeed, the photon has been suggested as an alternative 
to the candela as the definition of luminous intensity |f36l|. We use twin-beam states and their 
perfect photon-number correlations to measure the efficiency of a photon-number-resolving de- 
tector. This type of state can be produced at a wide range of wavelengths from many different 
types of source including nonlinear crystals, optical fibers, periodically poled waveguides, and 
atomic gases. This method has the advantage that it only assumes perfect photon-number cor- 
relation and does not assume the state has a specific photon-number distribution, nor even that 
it is pure. Despite these seemingly detrimental assumptions, the efficiency estimation presented 
has a large amount of redundancy leading to a relatively small absolute error of 0.4%. We show 
that this measurement is independent of the average photon number of the state, unlike the 
Klyshko method, making it more widely applicable. In particular, it is ideal for characterizing 
photon-number resolving detectors, and for use with bright reference states. PNR detectors are 



undergoing rapid development via a number of competing technologies. These detectors will 
play an important role in precision optical measurements and optical quantum information pro- 
tocols where photon-number resolution is necessary for large algorithms. Future development 
of direct efficiency calibration should focus on the issue of background, which can corrupt the 
state and thus the calibration. 
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